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The work described herein was done at the Research and Develop- 
ment Lahora-tories, Hn£;hes Aircraft Company, Culver Citj’", California, 
during the industrial terra of the Electronics Engineering curriculum. 
My control systera project v;as suggested by the Laboratories, and was 
carried out with their fullest cooperation as a part of a large scale 
study of airborne digital control systems. 

I shall assume that the reader is conversant with current termi- 



nology in the computer field (2). What few terms are not in general 
use I will explain. The words "order” and "instruction" are synonymous. 
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I. IIITROLUGTICII 



This thesis describes in detail a program and coding for a 
particular aircraft navigation display control system utilizing a 
digital computer, in less detail the equipment u.sed, and finally a 
number of suggestions for programming of control systems. 

Some modern control systems become veiy complex and require 
manj’’ individual computing or compensating elements. It seems likely' 
that a gain might be achieved by combining all these computing and 
compensating functions in one elaborate general-purpose computer, 
perhaps on a time sharing basis (l). Hither digital or analog tech- 
niques could be used. 

A digital computer offers certain inherent advantages in this 
application. Chief among these are, first, that changes in function 
of parts of the control system, such as adjusting compensation of a 
servo loop, are made simply as changes of numbers or instructions 
within the computer. Ho change of physical equipment is needed. 
Second, and more basic, is the memory, or storage of information 
that a digital computer contains. Thus, the control system can be 
influenced by its past performance, or by outside stimulation. Take, 
for example, an automatic factory. Tlien the output of the factory 
could be influenced by inventories and sales records which had been 
gathered and collated by itself. 

Such an elaborate control problem is presented in seme high-per- 
formance aircraft, and the computer on which the described program wa-s 
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run was a digital conputer for airborne use (3). 

Programoiing and coding a problem for a digital comouter have 
received much attention (4), but the approach has been one of re- 
ducing to a ninimim the human labor, either by admittedly inefficient 
computations (5,og.ll5) or by making the machine do the v/orl: of coding 
(6), This approach leads to overall efficiency v/hen a problem is done 
once or a times. If, however, a oroblem is repeated ad infinitum, 
as in a control system, human planning time starts to pay off as re- 
duced physical equipment due to the lessened storage capacity needed 
for the more efficient set of instructions. This v/ill be discussed 
later. 

In any machine that is very complicated some means of diagnostic 
testing should be provided. ‘This v;ill be elaborated.. 

Finally, a reasonably didactic descriotion of the program v;ill be 
given. This program was a first attempt, and it is on the basis of the 
mistakes made that the recommendations have been formulated. 
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The equipment used, while not the subject of this thesis, did 
influence the programming and coding, and so should be described. 
Peculiarities of the computer logic important to the program v/ill be 
described in later chapters when applicable. 

The computer used v/as a first model which had been used for some 
time for cii'cuit and component reliability and logical checks. Uo in- 
put-output equipment had been designed specifically for it, although 
input-output equipment v/as available that had been used on a very 
similar computer. One function of the Control Bos (Chap. II) was to 
adapt and interconnect the computer and the input-output equipment. 

The Brown Teleplotter was used intact except for changes in its in- 
ternal servo system to I'esoond to dc signals of the desired amplitude. 

The computer is of the general purpose, serial, binary digital 
type (3). It has an arithmetic unit, a control unit, and a magnetic 
drum memory unit. The a.rithmetic unit performs the arithmetic oper- 
ations of addition, subtraction, multiplcation and division. The 
unit consists of three one work circulating registers and a binary 
adder. The magnetic drum memory (3) provides one eight v/ord circulat- 
ing register, the "d" register, and fourteen sixty four v/ord storage 
bands. These bands are interpreted by the control unit as eight order 
bands containing instructions, and six number bands. One of the ntunber 
bands and one of the order bands are tied together, so that an instiuct- 
ion can be interpreted as a number and changed, written back into the 
memory, and then used as a nev/ instruction. Such changes are made as 
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the result of the computations in progress. The instructions avail- 
able to the coder are shown in Table II, 

The input-output equipment converts analog information (dc volt- 
ages) into an equivalent binary number at a sampling rate of once every 
fourteen magnetic drum revolutions, writes the most recently sampled 
binary number onto the drum, reads the most recentl;/’ vnritten output 
binary ntunber (v;ritten by the computer) and converts it to an equiva- 
lent analog quantity (dc voltage). It is multicriannel, having nine 
inputs and four outputs, v;ith one calibrating period ( 7 ). 

The plotting board consists of a chart and pen so arranged that 
two input dc voltages are interpreted and plotted as x and y deflect- 
ions of the pen on the chart. A potentiometer error detector feeds a 
modulator, whose ac output is proportional to the error. An amplifier 
then drives a motor in such a direction as to tend to make the error 22 ro. 
This is duplicated for the two axes. The speed of response of the pen 
is much faster than the contemplated plotting speed. 
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1. The Problem 

A navigation display vjas desired to demonstrate the ability of 
the Hughes Model I digital comouter to store much information for a 
long time. The operation of the device was visualized a,s follov/s. 
First, a map of the desired track is placed on a plotting board con- 
taining a movable clotting pen. The oen is placed in turn over the 
departure point and each of the enroute v/ay points, or "fixes". At 
each point a button is depressed which will enter the position into 
the memory of the computer. After a.ll points are entered, the pen is 
returned to the departure point and the computations started by throw- 
ing a sv/itch. The pen then "flys" to each of the fixes in turn at a 
fixed gro\ind speed on straight line courses. After getting to the 
last fix, or destination, the comouter vrould turn itself off after 
preparing the program for entering the next set of fixes. Provision 
was made for entering sixteen fix positions. A signal light would in- 
dicate v^hen the computer was ready to ha,ve fixes entered, or ready to 
fly. 

The Control Box had the follov^ing operational controls: tv;o hand- 

v/heels for Fast-Vest and llorth-South positioning of the pen during the 
entry portion of operation, a pushbutton for entering fixes and start- 
ing fly operation when ready, a function switch with Enter and Fly 
positions, and a switch to turn on recoi'ding pen on plotting board to 
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prevent smearing. Contained in the same hox were various inter- 
connections "between the computer, the injiut-output equipment, and 
the plotting board, various disconnect switches for the power and 
timing pulses, and reading and v;riting amplifiers required by the 
inpvit-output equipment. A schematic of this box is included as 
Figure 1 to indicate the functions, although the details are im- 
material to this thesis. 

It was desired to make the fly routine a feedback loop so that 
the present position of the pen could be sensed in the computations, 
and so that smoothing could be included by placing a lov; pass filter 
after the dc analog output. Smoothing could have been done in the 
program, and probably will be included at a later date by Eughes 
engineers. Smoothing in the program would insure, for example, that 
computational errors of a random nature wo\ild be smoothed out of the 
output, rather than being cumulative. Smoothing in the dc analog out- 
put v;ill accomplish the same result, but requires additional equipment. 

Smoothing should not be confused v/ith the essential filter required to 
\ 

smooth the sampling rate. 

In the course of the fly computations, the present position is 
written into the memory by the input equipment. The computer then 
reads this position and compares it with the fix position to which the 
**plane‘* is nov/ headed. By a trigonometric method involving the fixed 
speed towards the fix position, an incremental change in position for 
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the computing period is determined and added to the present position 
to give a future position. This future position is then fed to the 
output equipment, throxigh the low pass f ilter, and to the pen to posit 
ion it. This future position then becomes the present position for 
the next computation. The computations are repeated every 14 drum 
revolutions or 0.1 seconds, i.e., the sampling rate is ten cps. llote 
that if a computational error is made and the plane gets off course, 
it will still continue to fly to the fix on the nearest straight line 
course. IHien the plane arrives at the desired fix, the computer will 
change its instructions so that the next fix position will be used. 
After getting to the sixteenth fix position the computer v;ill stop, 
after setting in the proper instructions to enter nev; fixes. 

2. The Ikiuations 

The course and incremental position are computed by the following 
set of equations. 

Let x,y be the coordinates of the present position. 

be the coordinates of the next fix position, 
be the sampling period, 0.1 second, 
dx, dy be the coordinate distances travelled in tine dt. 

be the specified ground speed tov/ards the next fix 
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Then we can say that 

Yp -Y - tan 0 = dy/dx. 
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cos © = 1 / 1 / tan 

dx s Vdt cos 6. 



“S. 



dy I djc tan ©. 

/ cbc. 

Yn/ 1 = ^n / l^ese last two are the futuxe 

position, and become z, y for the next period. 

Difficulties arise from two sources, (l) Ambiguities due to 
sign of tan ©, and ( 2 ), the fact that the computer can handle only 
numbers v/hose magnitude is less than one. These v;ill be taken care of 
by using the conditional transfer clcS instruction. (Table II). ITote 
that no instruction is provided to specifically handle square roots, 
so an approximation will have to be made for cos 6. 

The equations finally used that follovf, resolve these difficulties 
See fig. 2. 

G-iven x, y, Xp, Vdt, Z. 

1. Pind ]Xp -■'■>■1. 

2. Jind |Yp -Yj. 

3. Yind lXp-Xl / iYp-Yl. Store for future use. 

4 . Yind |Zp-X|. - 1-Yp-Yi. Check sign. If 



5a. Positive. Yind 

(Yp-Y)/(Xp-X) = tan e. 

6a. Approximate cos 9 by polynomial. 

7a. Multiply Vdtcos6s dx. 



5b. negative. Yind 

(Sp-X)/(Yp-Y) c cot e. 

6b. Approximate cos 6 by poly- 
nomial . 

7b. Multiply VdtcosS = dz. 
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8a. Multiply dxtanQ = dy. | 

9. Divide dx/(Xp-X). Check sign. 
10a. Positive. Do nothing, 
dx z. dx*. 



8h. Divide dx/cot9 =. dy. 
If 

lOh. negative. Subtract 
0 - dx = dx*. 



11. Divide dy/(Yp-Y). Check sign. 
12a. Positive. Do nothing | 



If 

12b. Negative. Subtract 



I 0 - dy =. dy*. 

13. Add X ^ dx* s. . Write as output. 

14. Add Y / dy* - Yn+i •/rite as output. 

15. Subtract (\.2:p-Xl/\Yp-Y|) - z. Check sign. If 

16a. Positive. Start over l6b. Negative. Start over 

v/ith same Xp, Yp. with new Xp*l, Yp^l. 

17. Count up a number such that after executing 16b sixteen times, 
a change of sign occurs to stop computer. 

In the above, step 4 determines v/hether tan 6 is going to be 
greater than one. If so, cot 6 is used. Steps 9a3id 11 insure that 
dx and dy are of the proper sign. Step 15 checks how closely the 
ple.ne is to the fix position. V/hen this distance, or rather a dis- 
tance xirhich is slightly greater than the distance to the fix, becomes 
less than some small number z, the plane has arrived, and the next fix 
position is used. 

The polynomial used to auproximate cos 6 from tan 6 or cot 0 is 
derived by a Chebyshev polynomial operation, and is accurate to 0.3/’ 
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v;ithin the range used; 

cose i 0.14743:^-0.4340x^^0.9965, where x = tane, or 
cose - y(. 147^^-. 43407^ /0. 9965) , where y =. cot 9 . 

This approximation is easily done on the computer, and is easier than 
conventional algebraic methods would be. 

3. The Ply Program 

Using the instructions contained in Table II and these equations 
for a guide, the program was designed and is presented in Table IV. 
Table III shov/s a flow diagram of the instruction sequence. Before 
detailing the instructions some general explanations are required. 

On a machine v^hich has available only such unsophisticated in- 
structions as shown in Table II , the detailed coding must proceed 
closely with the program, since the order in v/hich operations are per- 
formed depends upon, for example, v/hen certain numbers become avail- 
able from the memory. It may be advantageous to delay some operations 
till a more convenient time. Thus the sequence of instructions may 
seem odd until the coding is studied in Chapter III. 

The subscript on instructions involving the D register denotes 
the sector number in this eight vrord memory. 

The B register acts as an accumulator for the arithmetic unit; 
tliat is, any number written into the B register adds algebraically 
to the number already there. In starting a new operation it is neces- 
sary to insure that the B register is cleared, as in instruction 4 
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Any instruction involving "froni B”, such as ha, hm, etc,, v/ill clear 
the B register. 

All other registers and the memory v;ill hold a number until a 
new number is written over it, replacing the old number with the new 
number. 

Instruction ifl will be called simply ,rl, etc. 

#1 clears the B register. 

transfers the fix Xp to the D register. T>As instruction vrill 
have to be changed as we orogress from fix to fix. if 3 does same for 
Yp, but v/ill be a series of instructions, as will be seen in Chapter III, 

7 ^ and ^ bring in the uncorrected values of present position, Xy 
and yy , from the memory. These numbers are then divided by a scale 
factor (7) to become the corrected x a,nd y, and are stored in D (up 
to #17). 

#18 to #38 are manipulations to obtain K'#L’, which is a rough 
measure of the distance to go to the fix, see Figure 2. Some inter- 
mediate results are stored in D for future use. 

The q.uantity K‘/ L' can be a number from 0 to #2. This computer, 
which is designed to handle numbers between -1 and /l (modulo 2), in- 
terprets a nmber beWeen 1 and 2 as a negative number. Thus chedcs 

this fact, and if negative, indicating that K'/L' is greater than 1, 
writes in its place an arbitrary positive number, 

#44 to jf49 form K'-L' as a measure of the octant of 0. If in an 
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octant such that tan 9 is less than 1, the coapxite.tions use tan 9. 

If tan 0 is greater than 1, the computations use cot 9. Shis avoids 
overflow on the follov/ing division. #50 senses this fact. 

#51 to #54 form tan 9, OR to #?6 form cot 9. 

^Vhichever form is used, #55 to #72 form either cos 0 or cos 9/cot 
9 by the approximation shown on page 10. 

#78 and #79 bring 1T*-L' into the A register again to check sign. 
This 'wa.s an easier operation to code than to store the actual decision, 
#50, previously made, i.e., the decision was made twice. 

#8l to #93 fora dx and dy for the case when cot 0 is used. The 
sign of dx and dy may be wrong. 

#94 to #102 form dy and dx for the case when tan 9 is used. Again, 
the sign may be wrong. 

#103 to #118 check that dx is in the same direction (has same sign) 
as X coordinate distance to the fix. If wrong, the sign is changed by 
subtracting from zero. 

#119 to #137 does same for d^-. llote again that a definite possi- 
bility for overflow exists in the tv/o divisions, v/hich are used to check 
parity of sign, when the present position approaches a fix. dx and dy 
must be larger than the actual change in present position due to the 
time lag effect of the smoothing filter. T}-'is mistake was corrected 
by utilizing a multiplication to check parity of sign in the actual 
machine, but v;as left in this progrsjn to illustrate that not only must 
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you think in the language of the computer, but also in terms of the 
overall machine considered as a control system, i.e., its analog as- 
pects. 

#133 to #149 form the future position and write this position 
into the memory in a position v;here the output conversion device can 
read it. 

#150 to 7fl5^ compare K’/L*, or roughly distance to fix, iirith a 
small arbitrary number, z. l/lien the pen has approached within z of 
the fix position ifl55 senses this change of sign. If not within z of 
fix, several waits are necessary so that recycling back to the ?fl 
start of the program will be after the l4th revolution of the magnetic 
drum. This is to insure synchronism with the sampling rate of the in- 
put-output equipment . 

If v/ithin z of fix, #l60 to #l65 take #2, add a number to the in- 
struction such that the "words to wait" (see Chapter III) are changed 
to give a new Zp and Yp, and write the new back in position so tliat 
it can be used the next time around. A^ach time a fix is changed a 
number is added to the count number such that after the sixteenth fix 
is passed the count number becomes greater than 1, overflows, and 
appears as a negative number to #175. irl79 then i/rites 0 in output 4, 
which is connected through the output device to a signal lamp, puttiTg 
out the lamp. Later we will see that the lamp had been put on after the 
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last fix v;as entered. ?^180 to #132 v;rite in #2 of Table ?, the 
entry program (see next section), so tha,t the instructions will 
be correct for its initiation. 

#183 uses the same order as some nunber not zero so that irl84 
will stop the computer after getting to the sixteenth fix, or the 
destination. 

4. The Sntiy Program 

In Table Y are the instructions used to enter fixes into the 
memory. 

'f2 to Tf6 change the "words to wait" in the same manner as in 
previous section. Since this operation is done before using the 
order, the first "old am instruction" will be some number, meaning- 
less as an instruction, v;hich v;hen added to v;ill be the desired "new 
am instruction". This v/ill be discussed further in Chapter III. 

if? to 7 fl 6 alter present position, v;hich in this case is the 
desired fix to be entered, by the same scale change as in previous 
section. 

#1? is an example of programming influenced by the coding re- 
quirements. It v;as necessary to enter this number into the same 
position in the memory as the "old am instriiction" had been. This 
required sixty four word tines, and so the intervening time might 
as v/ell be used for some more useful purpose. 
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#21 makes use of the number written by #17 as an instruction, 
and enters the fix position into the memory. 

Since the instruction follov;ing #22 must occur immediately a,fter 
writing the fix position (See Chapter III), and the fix positions are 
in separate places in the memory, #23 becomes a series of identical 
instructions, but \v*ith different "words to v/ait". Shis provides a 
means of determining v;hcn the sixteenth and last fix has been entered 
without utilizing a count number as was done in the previous section. 

The last of this series is called TfeB. 

Jor the first fifteen points, then, :/23 to #27 puts an arbitrary 
positive number in to 3, a zero into output 4 for the signa,! light, 
and stops the computer ready for the next fix to be entered. 

#28 to #35 v/rites the instruction which becomes it2. of the fly 
program in its proper place, vrrites a number v/hich is ne^’.rly full scale 
into output 4 to light signal light, v;rites a zero in the covint number 
used in the fly program, and stops the computer ready for flying. 

5. The Test Program 

The test program shovm in Table I v/as included for diagnostic 

testing of the computer components. It v/as of a tyoe used previously 

for testing the llodel I computer, and performs the arithmetic functions, 

various transfers, and utilizes all the number and order bands of the 

—7 —7 

memory. It does these operations on all binary numbers from 2~' to 1-2 '. 
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This program is not explained in detail since it is not the v/riters 
work, and is immaterial to this control system, ezcept in its oper- 
ational use. 
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Coding the prograni developed in Chapter II repnires a detailed 
knov;ledge of the comp-uter and hov; instractions follov/ one another. 

'Table I, contained in the envelope in the bach of the thesis, lists 
the coded problem. Programming and coding were done nearly simultane- 
ously. 

Three separate problems are coded in Table I : the fly program 

v/hose instractions are labelled 1 , the entry progrejn labelled and 
the test program labelled T. The instructions in each of the three 
are serially numbered v;henever possible. Table III shows a flov; dia- 
gram of instructions for the fly program. Reference should be made to 
Tables IV and V for the sense behind the coding. 

The Model I computer uses a relative address code (3) in v/hich 
the address of the ne:mt instruction is given by specifying the memory 
band in which the instruction occurs, and the nijmber of v/ord times tiH 
the nev; instruction is ready to be read. The system is relative because 
the sector address of an instruction (or number) is never given, ej:ceot 
for the first of a series of instructions, v/hich alv/ays is in sector 1 
in the selected memory band. (The m^^ sector is the m^'*^ word after the 
arbitrarily designated sector 64 in memory). 

An example v;ill follov/ to illustrate this point. Pirst , the "book- 
keeping" shorthand should be explained. The instruction itself is follov/- 
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ed by a series of three numbers in fable I. The first is the v;ords 
to v/ait after reading the order, and before executing the order. The 
second number specifies the number band of the memory involved if mean- 
ingful; thus, the instruction -f does not involve the memcr;’’, but am 
does. The third number specifies the memory band of the next instruct- 
ion, For examole, in sector 06, band 1, the instruction nd]_ 2,4,1 
3T occurs. Tliis means that the number in number band 4 should be trans- 
ferred to the D register, sector 1, after a 2 word v;ait, and that the 
next instruction, m^T, will be found immediately after execution in 
band 1. The instruction v/as then read out daring sector (word) ,/06, 
the computer did nothing during sectors ;/07 and ,/OS, tra.nsf erred the 
number found in number band 4, sector ^9, to the D register during 
sector v09, and road instruction 7r'4T during sector /lO. hote that the 
D regi 3 t^r is eight v/ords long, and arbitrarily sector 1 of the D 
register coincides with sector 1 of the sixty four word memory. Thus 
the D register repeats itself starting v;ith sectors 9. 17. etc. of the 
memory. 

The permissible words to v/ait before execution contained in any 
instruction can be anything from 0 to l6 words, except the wait in- 
struction (v/t) which can have 0 to 32 words wait. The v/ait instruct- 
ion is executed immediately, or rather it v/ould be better to say that 
execution corresponds to reading the next instruction. For example, 
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check instruction #1597 in sector 52, hand 2, which feeds into in- 
str-action ,fl7. 

ITote that a particular operation to he ijerforsied hy an instruct- 
ion nay he v/rong. As long as the address of the next instruction is 
correct the computer v/ill continue to cycle instructions blindly. 

This is of particular importance in trouble shooting. The usual pi'o- 
cedure is to make the instructions cycle, then worry about what mathe- 
matical results are obtained. 

The three programs all start in sector #1, but in different 
memory bands. The starting position is selected either by a set of 
buttons on the computer, or by the last instruction read by the con- 
outer. Thus, when instruction #l43y is read and the computer stops 
ready to start on the entry orogram, the address of the next instruct- 
ion is band 3* ^^hen the computer is started again by pushbutton, the 
instruction #12 v/ill be read first. 

It is by means of the variable words to wa,it portion of an in- 
struction that sixteen dif'^erent fix positions are entered and taken 
out of the nemoiy. Hov;evcr, in the acttial number representing a particu- 
lar instruction, the compliment of the v/ords to wait appears. Thus, 
the four digits representing the words to v;ait v/ill appear oddl^'' v/hen 
added to in the manner sho^vu. »’^’hen all four digits are zero, this 
means l6 words to wait; adding one least significant digit means 8 



( 19 ) 





• - 






i‘r^ 












4JP^s^sbsi, vi 



• 9^- 



I 



• » m, 

< vC 



# ^ I £^i# 

if- 







«•! SJ 



• I 



'^yi 



•« 



- - I 'r P 












f I 






!A| 



(m^ • » t ' 

. .^dlN4. «»•• .i 

* BMUM ^»«rM II ^ 0 





9l 



words to wait, etc. See ntunber band 2. 

Rather than taking up the detailed coding step by step, the 
interested reader may follow through as much of Table I as necessary 
to get the system in mind. Compare Tables I, III and IV x\rith emphasis 
on the operation of the conditional transfer instruction ckS and the 
variable v/ords to wait. 

It must be obvious to anyone who has followed the coding in any 
detail that there is much v/asted time and space here. The problem has 
become not one of electronics or control system design, but rather one 
for an ingenious bookkeeper and for nximerical analysis. Time did not 
permit completion. If this control system were tc be used in an im- 
portant application, hov/ever, such effort would be essential to a good 
design. 
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CMPi’SR 17. C-StTSEAL DISGUSSICIT A.TD IGCOIuCSLTIIlTIOlTS 



The test progran used in this device was a good one in that all 
functions of the computer were tested. However, diagnostic indication 
of troubles encoixntered was very limited. A test should give a positive 
indication to a serviceman of the location of troubles. The instructions 
available in even as unsophisticated a computer as the Model I should 
allow this by making more liberal use of the ckO instruction throughout 
the test program shown in Table I. Since a digital computer is so lit- 
eral and human coding so liable to error, one must know whether the ma- 
chine is at fault or not before spending the hours req.uired to trouble 
shoot a jprogram for a control system. 

Once assured of proper machine operation, several steps can be 
taken to get a program to cycle instructions correctly. As noted pre- 
viously, it is usual to cycle instructions, then run test problems to 
check mathematical operations. After entering a set of instructions 
into the computer in which a mistake is made, improper cycling will usu- 
ally result in circulation through an improper sequence of instructions. 

If this occurs, the following steps are of value with the Model I coraouter. 

1. Change an instruction, say about hialf way through the program, 
to a cld3 at some point where the computer vrill surely stop. If this 
half of the instructions cycles correctly and gives proper intermediate 
results from a test problem, put the instruction back and change another 
instruction say three quarters of the way through. Continue this pro- 
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cedure till the trouble is narrov;ed down. 



2. After narrov;ing dov;n, change one or more instructions to 
cycle only a sBiall portion of the orogram repeatedly, and get it to 
v/ork. This might be called reentrant checking. 

3 . Another method after narrov/ing down is to use test equipment 
built into computer (provided in most computers) v/hich v;ill allow exe- 
cution of a single instruction at a tine. C’leck each step in a suspect- 
ed region. 

4. In troublesome cases it may be necessary to check each in- 
struction with an oscilloscope. As might be suspected, this is time 
consuming and should be avoided if possible. 

5. A good deal of ingenuity helps more than set rales. 

bTien proper cycling is accomplished some previously prepared tost 
problems should be used to check operations. The method used v^ith this 
set of programs wa,s to disconnect the input equipment, write in known 
nimbers where the in out numbers xirould have been written, and check either 
the final output, or^the numbers held in the registers at some intejnedi- 
ate time, using the narrowing doTO technique previously described. 

If now the computations are properly made by the computer, the rest 
of the control system comes under scrutiny, especially the stability of 
the system and errors. Hote that the program in use is an integral -oart 
of the control system, and must be tho’ught of as part of the equipment. 

A useful concept in present control systems is the tiansfer function 
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either of the system or its comoonents. A prcgrara may he thought of 
in terms of a transfer f-unction involving principally functions of the 
term exp(-as), a time lag (9). Unfortunately, too little has been done 
to develop a theory for this concept. It is mostly a study of numerical 
analysis (9). 

lamping of a control system has an analog in a digital computer 
program. Derivatives and integrals are replaced in these discrete cal- 
c^alations v/ith differences and sums. lamping then becomes smoothing, 
v;hich may take the form of averaging or interpolation, or approximation 
of the required derivatives by differences. The problem is one of numer- 
ical analysis and will be discussed further only in connection v/ith noise. 

A simple method of estimating stability vras suggested by Dr. Jacobi 
for the control system described. Since it is a feedback system involv- 
ing positive feedback, oscillation could occur when the total phase shift 
became a multiple of 2pi. The fly program has a time delay of 0.1 second, 
or a phase shift of 01= U/lO, where w is the angular frequency. The RC 
filter used on the output equipment has a phase shift ^ 2 ~ wRC, v/here RC = 
0,5. Oscillation at w = 0 is meaningless. 0 - 0i / 02 ~ 2pi occurs 
when w s 15 pi, or at a frequency of 7.5 ops. Rov/ever, the RC filter 
has an attenuation of 2? db.at this frequency so no trouble was expected 
or encoiuitered. 1/hen the filter v/as disconnected for trial, oscillation, 
or at least random action did occur altho^lgh the pen v/ould of course not 
follov/ the ^rations of the numbers in the computer. ITote that the same 
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smoothing effect could have been attained v/^ithin the orogram had it 
been designed into the program. 

jSrrors v;ithin the computations (excluding a steady bias error) 
may arise from two principal sources. P.ound-off errors appear as 
v/hite noise. Large random errors may occur due to electrical inter- 
ference or errors in programming (see pg. 12), or due to marginally 
satisfactory eq.uipment within the computer. Kie random errors may be 
reduced by proper design and maintenance, but probably cannot be elimin- 
ated. Both types of noise can be reduced v/ithin the program by smooth- 
ing. White noise is most effectively reduced by some means of averag- 
ing or interpolation. Large random errors can be treated, like ignit- 
ion noise, by limiting the value of successive increments in an iter- 
ative process, e.g., limiting dx and dy in the fly problem. 

An error which occurs as a bias error, e.g., a wrong scale factor, 
can cause considerable trouble in a system involving feedback. Scale 
factor is contained in the input and output equipments, within the pro- 
gran, and in any analog components ejcternal to these, particularly the 
sensors. Thus the accuracy of the overall systen is still determined 
by the weakest link. To take advantage of the possible accuracy of 
most digital computers, then, attention should be directed to the de- 
sign of the analog components of the system, cs well as to the program. 

In a time-sharing system it would be of value to utilize one computing 
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cycle for recomputing and ;;riting a new scale factor into the program 
periodically, i.e., automatic calibration. 

Using a digital machine implies both quantizing and sampling; these 
should be consistent v;ith system design, ’i/hen high sampling rates are 
required, e.g. , high roll rates of present jet aircraft, a stringent limi- 
tation is placed on the program designer and the computer engineer to make 
an efficient (few v/ords per cycle) and rapid (many v;ords per second) ma- 
chine. Low quantizing error may be required in some applications. These 
tv;o conditions are tinfortunately contradictory in a serial machine. This 
is especially iinforttmate v;hen one equipment is used on a time-sliaring 
basis for many jobs. The use of a pai'allel machine may be required v/ith 
its higher speed of operation but greater amount of equipment, if the 
program designer is not dexterous enough. 

A program may be made more efficient by efforts of the numerical 
analyst in reducing the steps reqtiired to do the series of operations. 

The coding of a particular program is a related but different problem. 

In coding, each v/ait for the aopearance of a number in the menorj’’ is in- 
efficient, as are comoletely filled sectors when other sectors are rela- 
tively vacant (compare sectors 5 s-nd 55 of Table I). The requirements of 
coding will probably determine the details of the computer design in the 
long inn. Sfficient coding is the result of attention paid to a great 
many small details. 
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Words to 

S:/Tnl>oI Meaning: I^ecnte 



cn Transfer mmber in c register to specified 

menoiy band 1 

am Transfer mmber in a register to specified 

nTimber band 1 

dm. Transfer number in d register to specified 

number band 1 

bm Transfer number in b register to specified 

number band, b cleared. 1 

ma Transfer number in specified number band to 

a register 1 

me Transfer number in specified number band to 

c register 1 

md Transfer number in specified mimber band to 

d register 1 

ca Transfer mtmber in c register to a 1 

Aca Transfer number, disregarding sign, from c 

register to a 1 

da Transfer number in d register to a 1 

ba Transfer number in b register to a register. 

b cleared 1 

cd Transfer number in c register to d. 1 

ad Transfer number in a register to d. 1 

dc Transfer niimber in d register to c. 1 

bd Transfer nijmber in b register to d. 1 

b register cleared 

/ Add number in a to number in b, sum in b 

register 1 

- Subtract number in a from number in b, 

difference in b 1 

X Multiply a by c, result in b register l6 

div Divide b by a, result in c, b cleared l6 

v;t Do nothing for (1 to 32) vrords 0 

cfcO Check b register. If zero go to next order. 

If not, stop computer 1 

ckS Check sign of number in a.. If positive - - 1 

If negative - - 2 

Table II. Instructions available in Model I 
Digital Computer. 
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start Problem 



Check sisn of |Xp-Xl / Up-^l 



Check sign of jXp-x) - 



Check sign of IXp-Xl - lYp-S| 



-* 1 

I 

40 

I 

'43^ 

If 
i 

94 ^81 

/ i 

102v .93 

^103^ 

I 

^114v 

115 ^llb 

\ 

^119 

^lio. 

131 ^153 

i i 

132v .13‘f 

i 

^155 V Check arrival at fix 

15b ^160 Count up words to wait of order 

1 I Check serial number of fix 

-159w ^175v 

^176 ^177 Prepare for entry routine 

I 

184 Check zero. Computer stops 



Check direction of dx 



Check direction of dy 



Table III. Flow Diagram of Fly Instructions 
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moldcn^i 




X 




X-qX 


\ X-qX Jo Oslo ioeriD 




ec 


04^ 






Sl^ 






X-qX 


- X-qX Jo ftsXB sloeriO 


O' 


oe 


X5 






dT 




4^5 


t-qX 


- X-qX lo flsXa iloeriO 


X8 


08 








ce 


cox 


SOX 




x5 Jo aolSooilb jCooIO 


dXI 


wx 


5XX 






8XX 


exx 






’tX* lo aot^oeitb sloeiiO 


m 


OCX 


xcx 








ecx 


sex 


*io6*io Jo 
XU 

oaliuot 


X£J d-43 XBVX’I'IB 3tO©ri0 
oJ Bbiov qu SapoO 
Jo lecfniijfl XjsXiee :So©riO 

lol ©lBqe*i<3t TTX 


08X 

5TX 


52X 

8TX 


d5X 

eex ' 



BqoiB 'xed’uqiaoO .oies jloedO ■A8I 



enoiifowiitenl '^X'5 lo ma'isBia woX*3 .1X1 ©XcfaT 



8S 



Inst. 


i:o. in A xlP-c:. ,,o.in 2 Lsr. 'Jo. in C 2er. "o.in nt’ orv 


l.ba 


0 


2.na 




3.130 


'^P 


^.cd3 


in Do 


5. ad^ 

6. \\'t 


Xp in Ih 


7.md6 


^ii ip IK 


8.ma 


T-a 


9./ 




10. ma 


s.c. 


ll.div 


0 y 


12 . cdy 

13. d£a 


Y in By 


14./ 




15 -na 


S.C. 


l6.div 


0 X 


17 . cd^ 


X in D^ 


iS.d^a 




19./ 




20 . d^a 


X 


21.- 


Xp-.w r E 


22 .bd2 


0 K in D 2 


23.d2C 


K 


24,Aca 


II 


25 .ad]^ 


E‘ in Di 




Table IV. Progrejn of Ply Instructions. 
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Inst . 


llo.in A Re,°:. 


x.o*in B Her. 


26 . :ga 


Yp 




27. i 






28.dr,a 


V 

X 




29.- 




II 


30 . bd^ 




0 


31.d4c 






32.Aca 


L' =tLl 




33.ad5 






34.d^a 


!£' 




35. / 




K' 


36 . d 5 a 


L« 




37./ 




K'/L> 


33 . ba 




0 


39.ckS 


If/, to 4o 
If-, to 4l 




40 .v/t 






4l.na 


1 

4 




42. v/t 






43 .cn 






44.dia 


x» 




45 ./ 




Z' 


46 . d^a 


L' 




47 .- 




K«-L' 


48. ba 


K'-L' 


0 


49 . an 






50.ckS 


If /, to 51 





If to 73 



1 ^ 0 . in nemory 



L in Jk 



L* in 



K«/L' nr 4 

in n. 



Y • ' 

Xi. ^-U 



in m. 



I'able 17. Continued 
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K’o.in B i.ep- 



ITo.in G -le^. 



jTo. in menor^/' 



Inst . 


1*0. in A jxee. 


ITo.in B 


51 . Aca 


L' 




5Z.i 




L« 


53 • 


X 




54 . div 




0 


55 . ca 


tan©, cots 




56. adi 

57. x 




tan^ 


58 . ba 


tan^ 


0 


59. ad2 

60 . d2C 

61 . x 




tan^ 


62 . bd^ 




0 


63 . ma 


- 0.4340 




64.x 




-.434tan^ 


65.M4 




0 


66. d5c 

67 . ma 


0.14?4 




68.x 




.I4y4tan^ 


69.d4a 


-.434tan^ 




70 ./ 




sum 


71 . aa 


0.9965 




12.4 


next inst 77 


cos 0 


73.dia 


K’ 




74 ./ 




K« 


IS. Leo. 


L> 





L' /K=±.tan^ 



tan in D]_ 



tan^ in D 2 

tan^ 



tan^ in 



tan 



4 



- . 434tan^ 
in D4 



Ta’ole IV. Continued 
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Inst. 


llo.in A Eef". 


ITo.in 3 Re-p*. 


ITo. in e.Rejor. 


Jo. in rneinon^ 


76.div 


’Text inst 55 


0 


K'/L' = cot9 




77. bdi| 

78. wt 
79 


K»-L' 


0 




cos or cos/cot 
in D4 


SO.ckS 

8l.d4c 


if /, to 94 
if to 31 




cos/cot 




82 . d^a 


cot 








83.x 




cos 9 






84.bd4 




0 




cos in D4 


85. me 






Vdt 




86. dZ)^ 


cos Q 








87.x 




dx 






88. ba 


dx 


0 






89.ad2 








dx in I >2 


90.-/ 




dx 






91.dia 


cot 6 








92.div 




0 


dy 




93.cd5 


next ins 103 






dy in d^. 


94. me 






Vdt 




95.d4a 


cos Q 








96.x 




dx 






97. ba 


dx 


0 






93.ad2 








dx in I >2 


99. die 






tan 0 




100. X 




dy 








Table IV. 


Cont inued 
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Inst. 


ITo.in A Ree^ 


llo.in B Tier. 


101. M5 




0 


102. wt 






103. d3a 

104.i 


V 




105 . d^a 


X 




106.- 




Xp-X*K 


107 .'bd4 




0 


108. d£a 


dx 




109 . •/ 




(lx 


110. d^a 


K 




lll.div 




0 


112 . cd]_ 






113.dia 


dx/K 




114. ckS 


if /, to 115 
if to 116 




115 . wt 


next ins 119 




116 . d2a 


dx 




117 .- 




-dXi dx* 


113. bd2 




0 


119.d£;a 






120./ 




Yp 


121. d'^a 


Y 




122.- 




Yp-Y-L 


123. bd4 




0 


124. d^a 


dy 




125 ./ 




!3y 



ITo.in C Ee?. 






1^0, in TieTnory 



K in D4 



dx/K in D]_ 



dx* in I>2 



L in D 4 



Table IV. Continued 
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Inst . 


llo.in A Ref. 


llo.in B Ref. 1 


To, in C lief:, llo.in llemonf 


126 


L 






127 . div 




0 


dy/B 


128. cdi 






dy/L in I>i 


129 .d]^a 


dy/L 






130.clcS 


if /, to 131 
if to 133 






131.d2a 


<iy 






132 ./ 


nert ins 13 5 


dy 




133 • 


dy 






134 .- 




-dy=dy* 




135 . ba 


dy* 


0 




136. ad^ 

137 . / 




dy* 


dy* in 


133 . d'pia 


Y 






139 ./ 




Y/dy*--Yn/l. 




I40.bd5 




0 


in B5 


141. d2a 


ds* 






142./ 




dx* 




143. dga 


X 






144./ 




X/d-x* sXn/l 




145. ba 


Xn/1 


0 




l46 . d^c 






Yn/1 


147. am 






in output 


148 . adj^ 






in Di 


149 . cn 






Yny^i in output 


150 .ma 


K'/L' or i 







Table IV. Continued 
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Inst. 


bo. in A Re^r. 


llo.in B Iter.* I'o.in C Rer.» Ilo.in memory 


151.-/ 




K’/L* or 


152. ma 


z 




153.- 




(kVl’)-z 


154. ba 


(K»/L«)-z 


0 


155. ckS 


if 7^, to 156 
if to 160 




156. wt 






157. wt 






158. wt 






159. wt 


next inst 1 




160 ,wt 






161. ma 


old ma 
instruction 




162./ 




old ma 
instruction 


163. ma 


1/16 




164 . i 




new ma 
instruction 


165. bd4 




0 nev/ ma inst 

in D 4 


166. ma 


167./ 






168. ma 


coTxnt no. 




169./ 




new co-unt no. 


170. ba 


new count no 


0 


171. wt 






172 . d^c 




nev; na 
instruction 


173 •cm 




nev; na inst 
in meaory. 


174. am 




nev; count no 
in memory. 


175.cl:S 


if /,to 176 

if -,to 177 





lable IV. Continued 
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1 



1 ) 



1 



f 



last . jno.in A I:er~, I\o.ia B iier. ITo.in C Ss; 

176 . wt next insl59 

(recycles) 

177 . roa 0 

178 . wt 
179 -am 

iSO.ina new ara inst 

181. wt 

182. am 

183 . / new am inst 

184. clt0 stop computer 



Table IV. Continued 



llo.in memory 



0 in output 4 



new an inst 
in memory 



36 









-^<•1 

^ mm 

ll. • 
4f 




.4 i>r4^ 





Inst. 


llo *in A Eep*. 


I'To.in B 


-o.in C Her. 


To* in memory 


l.ba. 




0 






2. na 

3. / 

4.iaa 


old am 
instruction 

1/16 


old am 
instruction 






5. / 

6. bd/j, 

7. na 


Xu 


new am 

instruction 

0 




nev; am inst 
in 


S.mdg 








Yu ^8 


9.7^ 










10. ma 


s.c. 








ll.div 




0 


X 




12.dga 


Yu 








13 . ^^2 








X in D 2 


14./ 




^u 






15.na 


s.c. 








16. div 

17. d£a 
IS.dijji 

19 .v:t 
20 .wt 


X 


0 


Y 


new am inst 
in memory 



Table 7. Program of l^ntry Instructions. 
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Inst • 


iro.in A lie#-. 


jlo.in I! Rer. 


ITo.in C Rer:. 


•:Io,in /^emerv 


21. am 








X in nenory 


22.cn 








Y in memory 


23 . nc 






0 




24. ma 


1 

2 








25 ./ 




1 

2 






26 .cn 








0 in output 4 


27.ck0 


Computer 










stops 








23. me 






ma instr 




29 . ma 


1-2"'^ 








30 ./ 




1-2-7 






31 . cm 








ma inst in 


32 , me 






0 


memory 


33 . am 








full scale in 


34 . era 








output 4 
0 in count no , 


35.clcO 


Computer 









stops 



Table y. Continued 
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